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ABSTRACT 

We revisit the problem of the formation of dense protostellar cores due to ambipolar diffusion 
within magnetically supported molecular clouds, and derive an analytical expression for the 
core formation timescale. The resulting expression is similar to the canonical expression ~ 
T^/Tni ~ lOrg (where rg is the free-fall time and r^i is the neutral-ion collision time), except 
that it is multiplied by a numerical factor C(^co)! where //cO is the initial central mass-to- flux 
ratio normalized to the critical value for gravitational collapse. C(/ico) is typically -^ 1 in highly 
subcritical clouds (/Uco <SC 1), although certain conditions allow C(;Uco) 3> 1. For clouds that 
are not highly subcritical, C{fico) can be much less than unity, with C{iJ,co) — > for ^Uco — > 1, 
significantly reducing the time required to form a supercritical core. This, along with recent 
observations of clouds with mass-to-flux ratios close to the critical value, may reconcile the 
results of ambipolar diffusion models with statistical analyses of cores and YSO's which suggest 
an evolutionary timescale ~ 1 Myr for objects of mean density ~ lO^cm"^. We compare our 
analytical relation to the results of numerical simulations, and also discuss the effects of dust 
grains on the core formation timescale. 

Subject headings: diffusion — dust, extinction — ISM: abundances — ISM: clouds — ISM: 
magnetic fields — MHD — plasmas — stars: formation 



1. INTRODUCTION 

Ambipolar diffusion, the drift of neutral matter 
with respect to magnetic field and plasma, was put 
forth by Mestel &: Spitzer (1956) as a means by which 
magnetic flux can escape from interstellar clouds and 
thereby induce overall collapse and fragmentation. 
To explain the well-known inefficiency of star forma- 
tion, Mouschovias (1976, 1977, 1978) proposed in- 
stead that, rather than reduce the flux of a cloud 
as a whole, ambipolar diffusion gravitationally re- 
distributes matter within the central flux tubes of a 
magnetically supported cloud. He further suggested 
(Mouschovias 1979, 1982) that the formation of cores 
would occur on the ambipolar diffusion timescale 
r^j3 ~ Tg/Tni, where rg and t^i are respectively, the 
free-fall and neutral-ion collision times (discussion 
is also provided in Mouschovias 1987; McKee et al. 
1993; Galh & Shu 1993; and Mouschovias & Ciolek 
1999). For typical ion abundances within molec- 
ular clouds (e.g., Guelin, Langer, & Wilson 1982; 
Langer 1985; McKee 1989; Ciolek & Mouschovias 
1994, 1998), one finds r^^ ~ 10 rg over a wide range 
of densities. 

Numerical simulations in axisymmetric (Fiedler 



& Mouschovias 1993) and disk-like (Ciolek 
& Mouschovias 1994 [CM94], 1995; Basu & 
Mouschovias 1994 [BM94], 1995a, b) model clouds 
confirmed the theoretical scenario described above. 
No magnetic flux leaked out of these model clouds, 
while in the interior flux tubes supercritical cores 
formed through ambipolar diffusion; core-envelope 
separation ensued, in which dense cores collapsed dy- 
namically while embedded in massive, magnetically 
supported (subcritical) envelopes. For highly sub- 
critical model clouds (i.e., with initial central mass- 
to-flux ratio less than 1/e times the critical value), 
cores did indeed form (i.e., the central mass-to- flux 
ratio reached the critical value) on a timescale tcore — 
the initial central flux- loss time r^ ^^ = r^^ ^(,/2 (e.g., 
CM94, § 3; BM94, § 4). However, model clouds 
that had initial central mass-to- flux ratios {M/^b)cO 
that were less subcritical formed in a time tcore sig- 
nificantly less than r,j ^^ (e.g., model 5 of Basu & 
Mouschovias 1995b). This was also true for the 
recent model presented by Ciolek Sz Basu (2000; 
CBOO), which successfully reproduced observations 
of the mass and velocity structure within the LI 544 
prestellar core. For this model ^co = 0.8, where 
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;Uco = (M/$ij)co/(M/$B)crit is the initial central 
mass-to-flux ratio in units of the critical value for 
gravitational collapse, and tcore was found to be 



^ 0.2r^ 



This suggests that icore is dependent 



on the value of //cO- In this paper we derive an an- 
alytical relation for tcore(/^co) that is accurate for all 
subcritical clouds (i.e., ^co ^ !)• 

In addition to its general theoretical interest, this 
investigation is particularly timely, given that a re- 
cent compilation of magnetic field strength data 
(Crutcher 1999) shows that many clouds have mass- 
to-flux ratios close to the critical value, and some re- 
cent statistical studies of ages and lifetimes of proto- 
stellar cores and young stellar objects which seem to 
indicate that star formation can occur on a timescale 
~ 1 Myr (e.g., Lee & Myers 1999; Jijina, Myers, 
& Adams 1999) within objects of mean density ~ 
lO^cm"'^; that is, it occurs in ~ few rg rather than 
~ lOrg at these densities. Our general results in 
this paper (as well as the specific case presented in 
CBOO) show that these inferred lifetimes are compat- 
ible with the observed regions being slightly subcrit- 
ical and evolving due to ambipolar diffusion, in addi- 
tion to the possibility that they are already somewhat 
supercritical and evolving dynamically. 

2. DERIVATION OF THE CORE FORMATION 
TIMESCALE 

The fundamental equations we use were developed 
for disk-like model molecular clouds (using a cylin- 
drical polar coordinate system r, (j), z), and were pre- 
sented in § 3 of Ciolek & Mouschovias (1993; CM93) 
and § 2 of BM94; the results we derive here can be 
generalized to other geometries in a straightforward 
manner. Model clouds are isothermal, embedded in 
an external medium with constant pressure Pext > and 
contain neutral molecules and atoms (H2 with a cos- 
mic abundance of He), trace amounts of ions (atomic 
and molecular), electrons, and charged and neutral 
dust grains. The clouds are self-gravitating and mag- 
netically supported, with /ico ^ 1- 

As in Ciolek & Mouschovias (1996; CM96) we start 
by considering the motion of a neutral fluid element in 
a flux tube centered on the axis of symmetry (r = 0) 
containing mass M and magnetic flux ^b and inter- 
secting the equatorial plane at a Lagrangian radius 
Tj^j^t). In a frame moving with the neutrals, 






M d^B 



^% dt Jo 



M 



27r^-2-?'M-S^(?'M)^DK/)> 



:i) 



B 



where Bz is the magnetic field and v^ = Vi — Vn is the 
radial ion-neutral drift velocity; v^ = drj^Jdt is the 
neutral velocity and Vi is the ion velocity. In deriving 



this equation we have used the magnetic induction 
equation dBz/dt = —{l/r)d{rBzVi)/dr. The induc- 
tion equation reflects the freezing of magnetic flux in 
the ions, which is valid for the typical density regime 
in which cores form {n^ <^ 10^ cm~^; see § 3.1.2 of 
CM93). Core formation takes place in the innermost 
flux tubes of a cloud, and so we consider equation 
(HI) in the limit r^^ — s- 0. The numerical simulations 
cited in § 1 found that Bz was spatially uniform in 
this limit, therefore we may take ^b — '^Bzv'^ . In- 
serting this relation into equation (IT]), and dividing 
both sides by {M/^B)crit, yields 



1 djdc 
He dt 



2v^ 



1 



(2) 



Ambipolar diffusion in the central flux tubes results 
in Wp > 0, which means that plasma and magnetic 
flux are "left behind" as neutrals drift inward toward 
the symmetry axis (due to self-gravity, see below), 
and equation (||) shows that the dimensionless central 
mass-to-flux ratio fic will increase with an instanta- 
neous e- folding timescale r^, ^ = r^^/2t'j3(r^^). Given 
the definition of the ambipolar diffusion timescale 
T^P = r/v^ (e.g., Spitzer 1978, § 13. 3e; Mouschovias 
1979), it immediately follows that t^ ^ = t^^ ^/2. The 
behavior of r^ ^ during the evolution of various cloud 
models is depicted in figures 2d and 4d of CM94, and 
figures 2c and 6c of BM94; examination of these fig- 
ures indicates that r^ ^ normally decreases as a cloud 
evolves (i.e., as /Xc increases with time). 

A supercritical core exists for /ic > 1- Core forma- 
tion occurs at a time tcore when /ic(icore) = 1- Inte- 
grating equation (0) from t = to t = tcove gives 



_ [^ dfJ^ r 

Ecore — / , Tg>,c\l^c. 

Thus, we may write 



C(/ico)T< 



#,cO' 



where 



C{n. 



cOj 



McO /^C ''"*,cO 



(3) 



(4a) 



(4b) 



and r^j ^Q = r^, ^.(^co)- The factor C{fico) defined in 
equations (|4a| ) and (|4b|) is the initial flux constant of 
a cloud. In a sense, it is a measure of the "distance 
a cloud has to go" in its central mass-to-flux ratio 
to attain the critical state /Xc = 1, given a particular 
initial value /ico- The "closer" /ico is to unity, the 
smaller we expect C(/ico) to be. Since t^ ^{fJ-c) < t^ co 
for fic > Hco (as noted above), equation (^) provides 
the limiting relation 



C(/Uco) < In ( 



(5) 
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Thus, C(^co) ^ as /ico —" 1- 

To proceed further, we need an expression for 
T^^(/ic) to evaluate the integral in equations @ and 
(Eq). This can be obtained from the force equation 
(per unit area) for the neutrals, which is given by 
equations (28a)-(28c) and (50) of CM93. In the next 
subsection we first consider models that account only 
for neutral-ion collisions and neglect the friction due 
to grains. The effect of neutral-grain collisions are 
then discussed in the following subsection. 

2.1. Clouds With Only Neutral-Ion Collisions 

The numerical simulations referred to in § 1 gen- 
erally found that the evolution during the core for- 
mation epoch is quasistatic, with little or no accel- 
eration of the neutral particles. Ignoring the accel- 
eration of the neutrals, as well as thermal pressure 
stresses (which are also negligible at these densities 
and lengthscales) , the evolution of the neutrals is de- 
termined by a balance between gravitational and col- 
lisional forces: 



~ a^gr + 



CTn 



(6) 



where a^ is the mass column density of the neu- 
trals, gr is the radial gravitational field, and Tni is 
the neutral-ion collision (momentum exchange) time. 
Solving this equation, we find v^ = Tni\gr\- Whence, 
again in the limit Vj^j — > 0, 



1 



$,c 



21-0 ('^m) 2(|fifr|rn 





where Tgr = {'f/\gr\)^''^ is the gravitational contrac- 
tion timescale (Ciolek &: Konigl 1998), which is ~ rg. 
(Note: in both CM94 and BM94, Tgr is referred to as 
the dynamical timescale r^yn-) 
For Tgr we may write 



A, 



/gr 



(c™„„. 



gco 



1/2 



(8) 



where G is the gravitational constant, m^ is the mean 
mass of a neutral particle (= 2.33 amu for an II2 



numerical constant that is solely dependent on the 
assumed geometry of a model cloud. For most ge- 
ometries, Ageo ~ 1. 
The neutral-ion collision time is expressed as 



1.4 1 + 



m^ 



1 



rui J ni{aw)in' 



(9) 



where m^ and m-i are the masses of H2 molecules 
and ions, respectively, ni is the ion density, and 



{aw)in is the ion-neutral collisional rate, which is 
~ 1.7 X 10-9 cm^ s-i for H2-Na+ and H2-HCO+ col- 
lisions (McDaniel & Mason 1973). The factor 1.4 in 
equation (P) accounts for the fact that the inertial 
effect of helium is neglected in the neutral-ion colli- 
sional force; for further discussion on this point we 
refer the reader to § 2.1 of Mouschovias & Ciolek 
(1999). 

Finally, we parameterize the central ion density 
through the relation ni.c ex; n^^, where, in this anal- 



ysis, we assume that /c is a constant. 



mg. 



k 



k{n 



n.c;; 



with 0.5 < A; ^ 



Strictly speak- 
1 for densities 



n 



n,c 



^ 10^ cm"^, as shown in numerical simulations 
(CM94) and analytical studies of the effect of ambipo- 
lar diffusion on ion abundances in contracting cores 
(Ciolek & Mouschovias 1998). However, as we shall 
see below, our neglect of the dependence of k on the 
density leads to only a minor error in the calculation 
of the integral in equation (|4b| ). Combining this re- 
lation with equations ^ , (^)7and (^ , and using our 
parameterization for nix(«n,c) discussed above, yields 



i-fc 



"-n.cO 



Ur, 



(10) 



Our remaining task is to determine nn,c(/^c)- In 
the disk model clouds, the mass density /On(r) [= 
"T'n^n(^)] is related to the column density cr^(r) in 
a magnetic flux tube through balance of thermal- 
pressure forces and self-gravity along the z-direction 
(i.e., along the direction of the magnetic field): 



PnC — Pcxt + izGcr^ , 



(11) 



C being the isothermal speed of sound (for a deriva- 
tion, see § 3.1.1 of CM93). It then follows that 



'^n.cO 



nr, 



Pext + 1 



-Pext + (crn,c/crn^ 



cOJ 



(12) 



where Pext = '^Pevit/'^Ga'^^Q is the dimension- 
less external pressure, in units of the initial self- 
gravitational stress in the central flux tube. Now, in 
the limit r^^ — > 0, {M/^b)c = (^n,c/Bz,c, and, there- 



gas with a cosmic helium abundance), and Ageo is a ' 



A^c 



(M/$ 



B) 




(13) 

In the last equality in equation ([l^) we have used 
the fact that the central magnetic fleld strength B^^c 
changes very little during the core formation epoch 
(Fiedler & Mouschovias 1993; CM94; BM94). That 
is, during this stage of evolution, ambipolar diffusion 
allows magnetic fleld lines to remain essentially "fixed 
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in place" (Mouschovias 1978, 1979). Inserting equa- 
tion (p^ ) into (p^), and then substituting into ( p^ ) 
and (pq), we have 



C(/Uco) 




l-k 



(14) 



where y = /^c/AicO- For a given value of Pcxt, C(/ico) 
can be obtained by integration of equation (|14|). Note 
that for A; = 1, equation (|l^ ) readily yields 



C(/i, 



cOj 



In 



1 

Afco 



(15) 



regardless of the value of Pext • Similarly, for the case 
k = 1/2, the integral ( |l4|) has the solution 



C{n, 



cOj 



Pcxt + 1 
-Pr-xt, 



1/2 



In 



pl/2 



(Pext + l) 



1/2 



^ext /^cO + (-Pc; 



t/^cO + 1 



1/2 



(16) 



The previously published model clouds cited in § 1 
typically adopted a value Pcxt = 0.1, appropriate for 
an isolated, gravitationally bound cloud. For such 



models the integral (14) may be solved by taking the 
limit Pext -^ 0, yielding 



C{^l. 



1 



cOj 



,2(l-fc) 



2 (l-k) 



for k < 1 



(17) 



(for k = 1/2, this expression is also easily obtained 
from eq. |16| in the same limit), and for k = 1, 
C(/Uco) is given by (|l5|) above. Figure 1 shows C(^co) 
as a function of //cOj for models with Pgxt = 0.1 
and k = 0.5, 0.6, 0.75, and 1.0, respectively. The 



solid lines are the solutions given by equations (15) 
and (|l7|), which assume Pgxt = 0. The curves with 
crosses, for the models with k = 0.5, 0.6, and 0.75, 
represent the solution for Pcxt = 0.1, obtained from 
equation (|14|) or direct numerical integration. Ex- 
amination of this figure reveals that the analytical 
approximation (O) for these small Pext models is in 
excellent agreement with the exact integration. 

Also shown in Figure 1 are the results for C{nco) = 
tcore/T^^Q as fouud in some of our earlier published 
core formation and collapse numerical simulations: 
model A of CM94 (Pcxt = 0.1, HcO = 0.26, k ~ 0.65 — 
see Fig. 2c of CM94), model 2 of BM94 (Pcxt = 0.1, 
;Uco = 0.324, k = 0.5), model 5 of BM95b (Pcxt = 0.1, 
;Uco = 0.5, k = 0.5), and the L1544 model of CBOO 

''That clouds may have mass-to- flux ratios near the critical 
molecular clouds (e.g., Larson 1981; Myers 1983). For instance, 
in 14 different clouds and showed that the ubiquitous relation 
a cloud) could arise from turbulent or nonthermal linewidths 
strengths comparable to the critical value for collapse. 



(Pext = 0.1, /ico = 0.8, k ~ 0.7). Comparison of 
the values found in the numerical simulations with 
those derived from our expressions for C(/ico) (eqs. 
01 "1113) ^^®° s^°^ excellent agreement. 

Notice that, for highly subcritical clouds (//co ^ 1), 
C(/Uco) ~ 1 for most realistic values of k. Hence, for 
these clouds, icore ~ '^^co' ^ originally argued by 
Mouschovias (1979, 1982). Models with fe ~ 1 and 
;Uco *C 1 are an exception: for these models t^ ^q does 
not decrease with increasing density, since the degree 
of ionization Xi^c = ''^i,c/''^n,c on n^~^ in the central 
flux tubes remains essentially constant as the cloud 
evolves. Under these conditions, r^ ^^ is a true e- 
folding timescale. 

However, C{fj,co) ^ 1 as ;Uco approaches unity, for 
all values of k. This region is of particular interest, 
given the recent OH Zeeman observations of mag- 
netic field strengths in dense molecular cloud gas by 
Crutcher (1999). The estimated mass-to-flux ratios 
are close enough to the critical value that a clear ob- 
servational distinction between subcritical and super- 
critical cores cannot be made. However, when the 
likely effects of random orientations of cloud axes 
and magnetic field projections on the plane of the 
sky are accounted for, Crutcher's inferred values are 
consistent with 0.5 ^ He ^ 1-5 (e.g., see Table 1 of 
Shu et al. 1999). Although his sample undoubtedly 
contains some cores that have already become super- 
critical (and therefore are in the early stages of col- 
lapse) , these values may also be indicative of ^co for 
typical star- forming molecular clouds. In this case, 
clouds may often be subcritical (and magnetically 
supported), but only moderately so. ^ From equa- 
tions dH) - (0) it follows that, for (1 - /ico) < 1, 
C(/Uco) — 1 — fJ'cO- Hence, C{nco) ^ 1 will be expected 
to occur for clouds with /ico values in the range sug- 
gested by the Crutcher (1999) data. This means that 
clouds with these values of /^co ^^^^ evolve and form 
cores by ambipolar diffusion on a timescale much less 
than r^^^o . 

Let us evaluate the magnitude of this effect numer- 
ically. From equations (0), (P), and (P) we find 



7.0 X lO^'A 



geo 



Uo-7 



yr, 



(18) 



where we have normalized Xi^cO to a value typical 
for cosmic-ray ionized gas at density n ~ lO^cm"^. 
Since, according to our results displayed in Figure 
1, C{fico) ^ 0.5 for clouds with fico ^ 0.5, we con- 
clude from equations ( ^a] ) and ([l^) that tcorc ^ 4 

value is consistent with the well-known linewidth-size relation for 
Mouschovias & Psaltis (1995) analyzed the data for the linewidths 
A.V oc R^'^ (where Av is the linewidth and R the radial extent of 
resulting from hydromagnetic waves in clouds with magnetic field 
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Myr for moderately subcritical clouds. This effect 
had already been seen in the L1544 model of CBOO, 
which had /^co = 0.8, resulting in C(/ico) = 0.2, and 
icore = 1-3 Myr (see their § 3). Thus, recent statis- 
tical analyses suggesting that star formation occurs 
in some regions on timescales ~ 1 Myr may he ex- 
plained by ambipolar diffusion models and observa- 
tions of subcritical clouds with central mass-to-flux 
ratios slightly below the critical value for collapse. We 
also note that the much longer lifetimes of a few x 10 
yr for giant molecular clouds as a whole can be ac- 
counted for by their lower mean density, yielding a 
higher ionization fraction Xi, particularly when ion- 
ization due to background ultraviolet starlight is ac- 
counted for (e.g., Ciolek & Mouschovias 1995). The 
higher values of Xi can greatly extend the lifetime of 
lower density regions even if the mass-to-flux ratio is 
close to critical. 

In the limit Pgxt S> 1, corresponding to model 
clouds situated in a high-pressure external medium, 
appropriate for cloud fragments or sheets located in 
a dense or massive cluster or cloud complex, or in 
a particularly hot medium or intense radiation field 
(perhaps suitable for proplyds), the integral in equa- 
tion (|l3) is found in the limit Pext ^ oo to be the 
same as that given by equation ((l^). For this case, 
clouds with //co < 1/e have tcore > r^^^o since r.^.^^ 
is again an e-folding timescale. This occurs in this 
limit because the mass density pn given by equation 
(pT|) is insensitive to the amount of mass within a flux 
tube (given by the column density a^)- This is seen 
from equation (^), which implies that the fractional 
variation of the density (5nn,c/?T-n.c is related to the 
fractional variation in column density (5(Tn,c/'7n,c by 

Snn,c _ r,('7n,c/q"n,co) ('^'7n,c/o"n,co) 



nr, 



(Cn,c/Cn,co) + -fext 
(/Wc//ico)^ + ^ext 



(19) 



where the last approximate equality is due to 
equation (|l^). For Pgxt —^ 0, it follows 

from (|l9|) that (Jnn.c/^-n.c — > 25(Tn,c/<7n,c ~ 

25fic/lJ'c, i-e., the central density increases twice 
as fast as the column density and mass-to- 
flux ratio. In the opposite extreme, Pcxt -^ 

OO, (5nn,c/?T-n,c ^ 2(o-n,c/o-n,co)(5<7n,c/o"n,co)/-Pcxt ~ 

2(/^c/Aico)('5/^c/^co)/-Pext -^ 0. For these models, nn,c 
is determined primarily by vertical confinement (com- 
pression) of the disk by the external bounding pres- 
sure Pext; and, as a result, there is little change in the 
density as the column density and mass-to-flux ratio 
increases in the central flux tubes. Since the den- 
sity does not increase much in this case, neither does 
the degree of ionization and r^ ^ is effectively con- 



stant as the cloud evolves; this follows from equation 
(|1^), which yields St^_Jt^_^ = -(1 - A;)(5nnx/"-n,c- 
(As noted earlier, k = 1 also yields 6t^ ^ = 0.) 
Figure 2 displays C(^co) for two model clouds with 
Pcxt = 20; one has A; = 0.5 and the other k = 1. The 
high-pressure analytical result C(/Uco) = ln(l/^co)) 
is a good approximation for these models. In sum- 
mary, we have demonstrated that ambipolar diffu- 
sion driven evolution actually occurs less rapidly in 
regions of high external pressure, since the density, 
which is then determined to a great extent by the 
bounding pressure, increases less rapidly as the cen- 
tral mass-to-flux ratio increases. 

Finally, before ending this subsection we make an 
important physical distinction. As described above, 
our new analysis shows that, for clouds with /ico ~ Ij 
C(/Xco) = tcore/T<s>co ^ -*-• Despite the fact that a 
core forms at a faster rate under these conditions 
(when compared to the uncorrected initial flux-loss 
timescale) , it should not be considered to be dynami- 
cal. Recall that our results are obtained using a qua- 
sistatic approximation for the equation of motion of 
the neutrals during the core formation epoch, con- 
sistent with behavior seen in our detailed numerical 
simulations. Therefore, although the timescale for 
core formation in this situation is much less than that 
which would have been predicted by the "classical" 
expression for a core forming due to ambipolar dif- 
fusion (i.e., tcore — "^4, co)' i^ should not be taken to 
mean that the mass in the inner flux tubes is under- 
going dynamical or accelerated collapse during this 
phase of evolution. As described in detail in the nu- 
merical simulations cited in § 1, magnetically-diluted 
collapse occurs within a supercritical core only after 
it forms, leading to star formation in an additional 
few Tff. 



2.2. Inclusion of Grain Collisional Friction 

We now consider the effect of neutral-grain colli- 
sions. Here we restrict our analysis to standard mod- 
els with Pext ^ 1- Detailed discussion of the effects of 
grains (charged and neutral) on the evolution of self- 
consistent model clouds is provided in CM93, CM94, 
and CM95. For typical grain sizes in dense molecular 
clouds (radii a ^ 10^^ cm) the bulk of the grains have 
a charge — e, where e is the electronic charge. More- 
over, as shown in CM93-95, rapid inelastic capture 
of ions and electrons on grains transfers momentum 
between charged and neutral grains, which can also 
result in effective coupling of the neutral grains to the 
magnetic field. 

We again make use of the fact that the evolution 
due to ambipolar diffusion during the core formation 
epoch is well-approximated as being quasistatic. The 
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motion of the neutrals is therefore given by 



~ CTnC/r + 



-^D + 



-\v^- 



(20) 



'ng 



The latter term on the RHS of equation ( |20[) is the 
force (per unit area) from neutral-grain collisions. 
The quantity r^g is the neutral-grain collision time 
and Ag [= {vg — Vn)/v^, where Vg is the grain veloc- 
ity] is the total magnetic attachment parameter of the 
grains, originally introduced in § 3.1.2 of CM93. Ag 
expresses the coupling of charged dust grains to the 
magnetic field: grains that are attached to magnetic 



field lines will have v„ 



Vi. 



and, therefore, Ag = 1; 



in the opposite limit of unattached grains (due to col- 
lisions with neutrals), Vg = v^ and Ag = 0. 



Equation (M) yields 



where a is the radius of the grains, ks is the Boltz- 
mann constant, and T is the gas temperature. 

Further progress is made by using the analysis of 
CM96, who studied the effect of ambipolar diffusion 
on dust-to-gas ratio during the formation of cores. 
As CM96 described at length, ambipolar diffusion can 
reduce the abundance of grains in a core when the in- 
wardly drifting neutrals "leave behind" the magnetic 
field lines and the grains that are attached to them. 
In particular, they showed that during the core for- 
mation epoch, Ag^cdfJ-'c/ fJ-'c = —dxg^c/xg^c- Therefore, 
when Ag^c 7^ (corresponding to some degree of at- 
tachment of grains to the magnetic field), the relative 
abundance of grains decreases as the central mass-to- 
fiux ratio increases. ^ Using this relation, along with 
equation (24), the integral (p3|) becomes 



^D(^Af) 



\9r\Tn 



1 + (Tni/'rng)Ag 



As in the preceding section, in the limit r^^^ 
follows from equations (|2|) and ( pi|) that 



<1>,C 



'-m +iK:Ag 



(21) 
it 

(22) 



the evaluation of the first 
the contribu- 



In the derivation of C(/ico) 

term on the RHS of equation (22) 

tion to the fiux-loss timescale due to the ions alone 

(which would occur for Ag = or Tai/r^g ^1) — 

follows exactly as presented in § 2.1 above. We focus 

now on the second term on the RHS of equation ( p^ ) . 

To determine its effect on C(/Xco) involves evaluating 

the integral 



1 du' r^ 

McO /^c ''"ng.c 



g,c 






/^gr.cO 



'"ng,cO . 
''"ng,c 

(23) 



g,c- 



Simplification comes about since 

I '''gr>c \ ''"ngjCO ^^11,00 '^gjC 



^gr,cO 



'ng,c 



n 



n,c ^g,cO •^gjCO 



(24) 



where ng is the density of grains and Xg = Ugjn^ is 
the relative abundance of grains. In the reduction of 
equation (^) , we have used equation (^ and the fact 
that the neutral-grain collision time 



rrin + m 



g 



1 



'ng 



rria 



ng{aw)gn 



(25) 



nig (^ rrin) is the grain mass, and {aw)gn is the grain- 
neutral collision rate, equal to 'Ka^{^kBT /irm^Y''^ , 

^This particular result is unique to our ambipolar diffusion 
discussed in CM94 and Ciolek & Mouschovias (1998). 



1 dj^T^ 

ll' T ^^'"^ 

McO H-c 'ng,c 



■A, 



gr 



'""S/cO 



■^giCf 



(26) 



where Xg^cf is the abundance of grains in the central 
fiux tubes at time icore, /^c = 1- CM96 (eq. [11]) also 
showed that 



■^giCO 



/^cO 



l + Al 



1/4 



1 + Alfi* 



(27) 



cO, 



where ^g is a dimensionless quantity whose inverse is 
related to the strength of the attachment of charged 
grains to the magnetic field. Specifically, equations 
(7), (8), and (9a)-(b) of CM96 yield 



A - 1 

1 + ^g^cO 



and therefore, for ^„ 



0, A 



g,cO 



(28) 



1 (grains 



are frozen to magnetic field lines), while, conversely 



for Act 



oo, 



A„ 



,cO 



(grains are completely 



unattached) . Numerically, 



./H.g 



2.32 X 10 



-2f B 



z,cO 



V35 /iG 

rrin 
2.33 amu 



10-^ cm, 
1/2 /lOKXl/2 

V T 



(29) 



(CM96, eq. [10]). For most situations of interest 
(10-6 cm ^ a ^ 10"^ cm, B^ cO ~ 10 /xG, and T ~ 10 
K), ^g « 1. 

We now have the necessary information to derive 

C(/Uco) in models clouds including the effect of grain 

models. It also affects the ion chemistry in contracting cores, as 
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friction. Combining equations pq), (|2^), and (pq)- 
(p3), along with the results for the ion contributions 
described in § 2.1, we have, in the limit Pext '^ Ij 



C(^co) 



1 ,2(i-fe) / 



l-^^o T^^ 



1+-42 



?^cO 



1/4 






CO V+-4iM^o 



(30a) 
if A; < 1, 



V McO / \-rngJ 



cO 



, 1+^g 



MJo 



1/4 



^+(-nJeo(l+iiA'^0 



(30b) 
if A; = 1. 



Examination reveals that, in the limits (Tni/rng)co — > 
(dynamically unimportant dust grains), or Ag — > cxo 
(grains are not well-coupled to field lines and move 
with the neutrals, thereby not retarding contraction), 
both equations ( |30a| ) and ( |30bD revert to the ex- 
pressions found in § 2.1 due to the effects of ions 
alone, namely, equations ( p^ and (|l5|). Further- 
more, for most circumstances, because 10~^ cm ^ 
a ^ 10~^ cm, resulting in ^g <§; 1 (see eq. p9|), 
clouds with 1 — fico ^ 1 have C(/ico) ~ 1 — ^cO) just 
as we had found earlier. 

Our derivation can be compared to numerical sim- 
ulations of model clouds that account for the dynam- 
ical effects of grains. In model B of CM94 (see their 
§§ 3.1 and 3.3), ficO = 0.26, B^^^o = 35 ^G, T = 10 
K, Xi^co = 2.0 X 10^"^, Xg,co = '7.1 X 10"^\ and the 
radius of the grains a = 3.75 x 10^^ cm. Insert- 
ing these values into the expressions for rni,cO (eq- 



)- 



'ng,cO 



(eq. 



) , and Ag (eq. [^ ) , and tak - 



ing k = 0.64 (see Fig. 4c of CM94), equation ( P0a| ) 
gives Cifico) = 0.86. In actuality, CM94 model B had 
C{fj,co) = tcorc/T^co = 0.92. Hence, our analytical ex- 
pression for C{fico) is again in very good agreement 
with the results of detailed evolutionary calculations. 
The core formation timescale tcore = C{fj,co)T^^o- 
Numerically, when the collisional effects of grains are 
accounted for, t^ ^^ is again given by equation (|iT 
but multiplied by a factor 



1 + 



Tni 



1 



^ ^'^S / cO •^g/^cO 



1 + 2.6 




where equations (^ and (25) were used to evaluate 
Tjii^cO and rng,cO) and we also made the approximation 
Agfi'^Q <C 1. In equation (^l|) we have used nig = 
Airpsa'^/S, where ps is the density of the solid grain 
material, and also the relation Xg^cO = ("T-n/"T'g)Xg! 
where Xg is the dust-to-gas mass ratio normalized 
to the canonical value for the interstellar medium. 
From equations ( |30a|) and (|30b|) it follows that clouds 
with /ico > 0.5 still have C{pco) <^ 1, and therefore, 
icore < T^,co- ^^T the values of r<j,^(, given by equa- 
tions ( p^ ) and (pl|), depending on the properties of 
grains in those clouds, values of tcore in the range 1 
to 5 Myr are still likely. 



3. RELEVANCE TO OBSERVATIONS 

Our results reveal that dense cores may be regarded 
as an ensemble of objects which are closely clustered 
on either side of the subcritical-supercritical bound- 
ary, as implied by recent observations (Crutcher 
1999). The subcritical objects are still evolving due 
to ambipolar diffusion, and the supercritical objects 
have crossed the critical barrier and may show ex- 
tended infall (e.g., the case of L1544; CBOO). This 
picture seems consistent with the observation that 
only some cores show extended infall, and that even 
LI 544 requires significant magnetic forces to keep the 
infall speeds to subsonic levels (CBOO), given that it 
has little turbulent support. However, the mean life- 
time for all the cores, from statistical studies, can still 
be ^ few Myr, due to the relatively short duration of 
the subcritical phase at the observed densities. 

Statistical studies of dense cores, usually defined 
as regions in which the mean density n^ ^ 10^ cm~'^, 
since they excite the NH3 (J, K) = (1, 1) line, reveal a 
lifetime ~ 1 Myr for starless cores (Lee & Myers 1998; 
Jijina et al. 1999). There are considerable uncertain- 
ties (at least of order one) in this determination, since 
the lifetimes are based upon a bootstrapping from 
inferred pre-main-sequence (PMS) lifetimes; the rela- 
tive statistics of cores with embedded infrared sources 
to PMS sources yields a lifetime for the embedded 
phase, and the relative statistics of cores with and 
without embedded sources yields a lifetime for the 
starless phase. This is done under the assumption 
that the formation times of the cores are uncorre- 
lated and uniformly spread in time up to the present. 
There is also a question as to whether the samples are 
complete, as many of them are biased toward cores 
that already contain stars, and regions that have al- 
ready undergone considerable star formation activity. 

A cautionary note regarding lifetimes comes 
from the submillimeter continuum studies of Ward- 
Thompson et al. (1994, 1999). In all starless cores 
detected in the submillimeter regime, the inferred 
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central density is in the range 10^ crn"^ — 10^ cm"'^, 
meaning that the cores have a much shorter gravi- 
tational contraction timescale than cores whose in- 
ferred central density is only 10'^ cm~^. In particular, 
at a central density nn,c = 3 x 10^ cm~'^, the spher- 
ical free-fall time r^^c = (37r/32Gpn,c)^/^ = 6 x 10"^ 
yr, compared to 3 x 10^ yr at a density lO^cm"'^. 
This means that if the submillimeter cores have a 
lifetime ~ 1 Myr, then they live for ~ 20 rg, making 
ambipolar diffusion very necessary at this stage (see 
Ward-Thompson et al. 1999, who even conclude that 
some ambipolar diffusion models evolve too rapidly). 
The resolution to this apparent discrepancy can only 
come if it is demonstrated that the submillimeter 
cores represent a small subset of the ammonia cores, 
and high resolution studies establish that most of the 
ammonia cores actually have central densities close to 
lO^cm"^. A complete sample of starless cores (with 
high angular resolution) is thus very necessary before 
drawing definitive conclusions about lifetimes. 

Another caution against preliminary conclusions 
about the timescale for protostellar evolution comes 
from a recent study by Kontinen et al. (2000), 
who measured the abundances of various molecular 
species in two dense cores: the starless, prestellar core 
Cr A C located in Coronae Australis, and the Cha 
MMSl core in the Chamaleon I dark cloud, which 
contains a class source. Comparing the results of 
their observational study with time dependent chemi- 
cal evolution models, they argue that Cr A C exhibits 
'later time chemistry' than that which is seen in the 
Cha MMSl core. Kontinen et al. suggest, as one pos- 
sibility, that the reason the Cr A C core appears to 
be chemically or chronologically older than the more 
physically evolved Cha MMSl core is due to the Coro- 
nae Australis cloud having had a longer ambipolar 
diffusion timescale than in Chamaleon I. Of course, 
based on our discussion in this paper, differences in 
the ambipolar diffusion and core formation timescales 
between these two parent clouds may reflect differ- 
ences in the values of //cO (and, thus, C(//co))) as well 
as disparities in xi^co and the initial central dust abun- 
dance Xg.co (which could result from variations in the 
local cosmic-ray ionization rate CcR and/or dust grain 
properties — see CM98 for a discussion on how the 
fractional abundances of charged particles depend on 
these particular quantities) , which would fix the value 
of r^ ^Q (see eqs. |l^ and |^) in each cloud. 

Finally, it is of interest to note that the recent Zee- 
man measurements cited above seem to indicate that 
the mass-to-flux ratio in some molecular clouds is 
comparable to the critical value for collapse. That 



it is so close to the critical value likely indicates the 
important role of magnetic fields in the formation, 
support, and evolution of molecular clouds. Unfortu- 
nately, since a rigorous and comprehensive theory for 
the formation of molecular clouds from out of the dif- 
fuse interstellar medium is still lacking (clearly an as- 
trophysical problem of some importance, yet outside 
the purview of our present study), ^ more definitive 
statements about the specific nature of /ico and why 
it takes on values ~ 1 for observed clouds cannot be 
made at this time. However, our lack of knowledge as 
to why yUco may take on certain values does not affect 
our analysis in this paper (nor our previously pub- 
lished models and parameter studies, cited in § 1), as 
we have formulated our result in such a way that it 
can be applied to all subcritical (/ico ^ 1) clouds or 
cloud fragments. 

4. SUMMARY 

In this paper we have reexamined the formation of 
protostellar cores by ambipolar diffusion in magneti- 
cally supported interstellar molecular clouds. Start- 
ing from basic principles, and utilizing the results of 
detailed numerical simulations, we have derived an 
analytical relation for the timescale for the formation 
of supercritical cores tcore within initially subcritical 
clouds as a function of /Xco, the initial central mass- 
to-fiux ratio in units of the critical value for gravi- 
tational collapse. The resulting expression is given 
simply by tcore = C(/Xco)r,j, ^g, where r^^^ is the ini- 
tial central flux-loss timescale, and the dimensionless 
factor C(/ico) is the initial flux constant of a cloud. 
Our analytical expressions are found to be in excel- 
lent agreement with the results of detailed numerical 
simulations of core formation by ambipolar diffusion 
in disk-like model molecular clouds, including a re- 
cent model presented by Ciolek & Basu (2000) that 
accurately reproduces the density and velocity struc- 
ture within the L1544 protostellar core. Our analysis 
has also included the dynamical effects of interstellar 
dust grains, which is also found to agree well with our 
earlier numerical models. 

We have derived simple, yet accurate expressions 
for C{iJ,co), valid for all ^Uco < 1. A strict upper 
limit C(//co) < ln(l//ico) is obtained. Clouds with 
;UcO -^ 1 generally have C(//co) ~ 1; and therefore 
tcorc ~ T^ co ' i'^ agreement with the early analysis by 
Mouschovias (1979). However, clouds with /Xco ^ 0.5 
have C(/Uco) *^ 1) thereby substantially reducing icore 
with respect to r^ ^^ . Clouds with these values of 
fico, which seem to be indicated in certain regions 
by recent OH Zeeman observations, have numerical 
values of tcore ^ few Myr. Hence, the evidence for 



^Interestingly, ambipolar diffusion can provide a means by which /ico can be deduced for molecular clouds. As discussed in CM96, 
reduction of the dust-to-gas ratio by ambipolar diffusion (see § 2.2, and eq. M] above) can be used to infer /ico from comparative 
observations of the relative abundances of grains in the cores and envelopes oTmolecular clouds. 
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comparably short timescales for protostellar collapse 
and evolution (e.g., Jijina et al. 1999; Lee & My- 
ers 1999) is not in disagreement with a picture in 
which observed dense cores are either somewhat sub- 
critical (and evolving towards the critical value due 
to ambipolar diffusion) or are already somewhat su- 
percritical and have started accelerated or dynamical 
collapse. 
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Fig. 1. — Initial flux constant C(/ico) = icorc/T^ ,.o for several model clouds with k = 0.5, 0.6, 0.75, and 1.0, respectively. Crosses 
represent the exact solution of equation ( p^ for Pcxt = 0.1 and k = 0.5 (eq. |lq|), 0.6, and 0.75. Solid lines are the solutions given 
by equations ( |l5[ ) or (Hj) for Pext ~ 0. Also shown are the values found from previously published numerical simulations: model 
A of CM94 (located by the center of ®), model 2 of BM94 (e), model 5 of BM95b (•), and the L1544 model cloud of CBOO (o). 
For most models with /j,cO ^ 1, C(/Xco) ~ 1- Above /Xco ^ 0.5, however, C{fico) is significantly less than unity. For 1 — /ico ^ 1, 

C{fJ.co) ~ 1 — ^cO- 



Timescale for Protostellar Core Formation in Magnetic Clouds 



11 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 T" 




Fig. 2. — Same as in Figure 1, except in the limit of high pressure, Pcxt ~ 20, for models with k — 1.0 and 0.5, given by 
equations, (Ha) and (M), respectively. 



